Nonequilibrium corrections in the pressure tensor due to an energy flux 
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The form P = a(u, v)U + b(u, v)JJ for the pressure tensor for a system submitted to an energy 
flux J, widely used for anisotropic radiation and proposed to be more general in |jj has been 
recently questioned by Nettleton 0. We provide a physical basis, in a completely different way, for 
this expression for anisotropic radiation and ultrarelativistic gases and we criticize some previous 
physical interpretations. We recall the necessity of an understanding of this kind of expressions in 
a thermodynamic framework. 



Radiation hydrodynamics P] is a subject of great interest in astrophysics, cosmology and plasma physics. However, 
the numerical methods proposed to solve the transfer equation for the specific radiation intensity X are in many cases 
computationally too expensive. Therefore, one usually considers the equations for the moments of X up to a given 
. order m ^|,|5| . Due to the dependence of the equation for the moment m on the moment m + 1 one needs to introduce 
a closure relation. If only the energy density u (m = 0) and the energy density flux J (m = 1) are considered, one 
must introduce a closure relation for the pressure tensor P (m = 2). The usual procedure is to introduce the so-called 
' Eddington factor \ defined by: 
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where U is the identity matrix, n := 4 and / the normalized energy flux defined as / := — . In the limit of isotropic 
radiation (Eddington limit), \{f = 0) = 1/3, while in the free streaming case x(f = 1) = 1. A number of different 
expressions for the Eddington factor have been introduced in the literature by interpolating between this limiting 
cases. Some of them have been obtained from maximum entropy principles. For instance, in J5|,^| radiation under an 
energy flux is studied by exploiting the entropy inequality, i.e. by maximizing a generalized flux-dependent entropy 
O ■ under a set of constraints, and an Eddington factor given by 



X= 3-3^4^3^ (2) 

is obtained. The same result is recovered in Q] from a information theoretical formalism, whereas different versions 
of this formalism have been used in jrUs] to obtain other variable Eddington factors. 
Thus, we observe that a flux dependent pressure tensor with the form 

P = a{u, J 2 )U + b(u, J 2 ) J J, (3) 

has been widely employed in radiative transfer. In order to obtain such a dependence, which apparently departs from 
local equilibrium, some authors jj],^]|8| have consider a flux-dependent generalized entropy and Gibbs relation. In 
addition, (|J) also appears in the study of an ultrarelativistic ideal gas under an energy flux by means of Information 
Theory @]. 

However, in spite of what is claimed by the authors in fl]||J|, (||) can be obtained for the two simple cases of 
radiation and an ultrarelativistic gas without abandoning the local equilibrium hypothesis. The equations of state 
and entropies appearing in (l],|^|| may be recovered as well. 
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First of all, let us notice that these systems are submitted to an energy flux, and not to a heat flux, because the 
condition of null global velocity has not been imposed. Therefore, it is not difficult to show that the considered 
situation corresponds to equilibrium (i.e., a purely advective energy flux), in contradiction to what is assumed in 
In fact, due to the symmetry of the energy momentum tensor of a relativistic system (i.e. = T Vfl ), the 
energy flux J verifies 

J = c 2 P. (4) 



This property can also be obtained directly for a system of ideal relativistic particles, with energy e, = \J m 2 c^ + pfc 2 
and velocity vl — c . Q p * . The total energy flux is 

^/pf+m 2 c 2 

J = ^2e l v l =:^2j i , (5) 
and introducing the expresions for e$ and v~l, it can be readily verified that 

ji = eivl = c 2 fi. (6) 

Therefore, equation holds for this system. 

Following (Q), the equations of state for a system submitted to an energy flux J (without any additional restriction 
on the particle flow) must be the same as the equilibrium equations of state of a moving system with momentum P, 
which can be obtained by simply performing a Lorentz boost to an equilibrium system at rest. Therefore, the systems 
considered in are nothing else but moving equilibrium systems. The distribution functions of both cases can 

be obtained by the use of Lorentz transformations as follows. The rest frame {Kg) equilibrium distribution function 
can be written as: 

f = e °o+/3o«o +tI ' ^ 

where g is the degeneracy, ao — — /3oMo and a = — 1 for bosons, a = 1 for fcrmions, a = for particles obeying 
Boltzmann's statistics and a = — 1,/xo = for photons and phonons. We consider the cases of radiation and a 
classical ideal ultrarelativistic gas, so eo = Poc. An observer at rest in a frame K moving with momentum — P and 
velocity — V with respect to the Kq frame measures an energy e = pc for a particle with momentum p (and velocity 
5) that verifies: 

eo = 7 (e " Vp) . (8) 

Substitution of (||) in (0) gives 

(9) 



oOc+f3e+Ipc 2 _|_ a 



where we have defined (3 := 7/3o and I :— —(3V/c 2 . Note that ((3, Ic) is the so-called coldness 4-vector. If we now use 
@ we obtain 

f e a+f3e+Ij +a ( 10 ) 

Now, we can recover the distribution function used in (IJ] for radiation by simply setting a = —1, a = 0, e = pc, j = ec, 
.9 = 2: 

whereas for the classical ultrarelativistic gas we obtain the distribution function proposed in (!]] setting a — 0. Once 
this distribution function is fully justified, the whole procedure in [Q holds. Thus, the results obtained in |l],|5|,[| are 
recovered and, in particular, defining the pressure tensor as the mean value of the operator 

N 
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and using (ffl), (0) is obtained. However, the physical interpretation given by this derivation is completely different to 
that given in fjjg||§. 

Clearly, the distribution of particles is anisotropic in the frame K due to the additional vectorial constraint (the 
global momentum P). The distribution function ( [j~l| ) allows the study of anisotropic equilibrium radiation (being 
the anisotropy due to the relative motion) but not the study of nonequilibrium situations. We propose the following 
heuristic argument to understand the physical situation. Although in (|J||| different methods were used to arrive 
to the equations of state of nonequilibrium radiation submitted to an energy flux, the authors never imposed the 
constraint of no global motion of the system. Therefore, when they made use of the condition of maximum entropy, 
they found an equilibrium moving system because equilibrium situations have the maximum entropy and the moving 
system verifies the imposed constraint of non zero energy flux. 

Let us remark another interesting feature of (|ll]) related to the physical meaning of temperature in this moving 
system. The distribution function (nl]) can also be viewed as a Planck distribution with an effective /3 e / given by 

/3 ef : =p + Ic cos e = P{l- —cos 0). (13) 

This expression is used, for instance, in cosmology in the study of the Cosmic Microwave Backgroung Radiation 
(CMBR) in order to take into account the relative movement between the Earth and the reference frame defined by 
the CMBR. By averaging over the angular dependence with the distribution function ( |TT| ) , it is obtained that 



so it is possible to define an effective mean temperature given by 

Te f := TIT = 7To ' ( 15 ) 

KBPef 

where To := k \ Q ■ Therefore, T e f is found to simply be the Lorentz transformation of To, according to Ott's trans- 
formation law M. This gives a simple interpretation for Ott's temperature, whose physical bases were controverted 
along the sixties . 

In Jj"l| , it has been argued that in some situations it is not possible to apply the methods of nonequilibrium 
thermodynamics for radiation. This is the case, for example, for two planar surfaces fixed at different temperatures 
Ti and Ti which exchange energy through a radiation field. Photons travelling in one direction are characterized by Ti, 
and the ones travelling in the opposite direction by T2, so it is not possible to assign a single temperature to radiation. 
Based on these arguments, recently ]2j has criticized the thermodynamical methods used in |||(| in their analysis of 
anisotropic radiation. Hovewer, we have proved that the system considered in |^|J^] is, in fact, an equilibrium moving 
system and therefore the criticism does not hold in this case. Let us remark that in the system considered by Essex 
jnj the distribution function is characterized by a double peak, whose relative heights depend on direction, while 
ours has a single peak whose position varies with direction. Thus, it is possible to define an angle-dependent effective 
temperature. 

In addition, in Q the validity of expressions like (||) for the pressure tensor has been questioned, both for gases and 
for radiation. We have seen that in the cases of equilibrium moving radiation or an ultrarelativistic gas, the pressure 
tensor adopts an anisotropic form due to the presence of an additional vectorial constraint (i.e. J). We think that 
these simple problems can serve as a guide to more complicated nonequilibrium situations. Therefore, it seems a 
plausible possibility that for a nonequilibrium system submitted to an energy flux and zero mass flow, the pressure 
tensor also depend on the energy flux, like in equilibrium. If that is the case, the dependence must have the form 
in ([l]) because, from a purely algebraic point of view, the most general tensor that may be built up in presence of 
a vector J must have the form a(J 2 )U + b(J 2 )JJ and according to the definition for the pressure tensor trP = u. 
However, this question remains open, and such a form for the pressure tensor is not free of difficulties, as pointed out 
in ||. Taking into account these criticisms, and the fact that expressions of the form (||) are widely used in radiation 
hydrodynamics, the convenience of finding a consistent thermodynamic scenario for these systems arises. We should 
also note that some variable Eddington factors \ have been proposed |?j using maximum entropy principles without 
a careful interpretation of the generalized flux-dependent entropies that naturally appear in the formalism. 

A plausible framework to understand these nonequilibrium flux-dependent entropies appearing in radiation transfer 
may be Extended Irreversible Thermodynamics (EIT) fl2|| . According to EIT, both temperature and thermody- 
namic pressure should be modified in nonequilibrium situations, if a generalized flux-dependent entropy function is 
considered. Up to second order in the fluxes, one has 
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s(u, v, J) — s eq {u, v) + a(u, v)J ■ J (16) 

and if pressure and temperature are defined, as usual, by the derivatives of the entropy function, one can easily obtain 
flux-dependent equations of state: 

§-? < 18 > 

where T is the kinetic or local-equilibrium temperature and p the local-equilibrium pressure and 8 and tt their 
generalized flux-dependent counterparts. In addition, the resulting pressure tensor was supposed in |lj to adopt the 
form: 

P = ttU + tpJJ, (19) 
where tp is determined by the requirement that trP = u. 
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